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Q-T Abstract 
CD 

We evaluate a Laurent expansion in dimensional regularization parameter 
e = (4 — d) /2 of all the master integrals for four-loop massless propagators 
up to transcendentality weight twelve, using a recently developed method of 
CNj ; one of the present coauthors (R.L.) and extending thereby results by Baikov 

and Chetyrkin obtained at transcendentality weight seven. We observe only 
multiple zeta values in our results. Therefore, we conclude that all the four- 
qq ! loop massless propagator integrals, with any integer powers of numerators 

and propagators, have only multiple zeta values in their epsilon expansions 
up to transcendentality weight twelve. 
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About one year ago Baikov and Chetyrkin published their results of eval- 
uation of all the master integrals for four-loop massless propagators in a 
Laurent expansion in dimensional regularization parameter e = (4 — d)/2 up 
to transcendentality weight seven 

These integrals are associated with graphs depicted in Fig. 1. The corre- 
sponding coefficients at powers of e turned out to be linear combinations of 
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(see, e.g., Ref. 0|). In our recent publication [3], we studied three four- 
loop non-planar massless propagator integrals (corresponding to non-planar 
graphs) where, according to the analysis of Brown [4J, one could meet not 
only MZV but also Goncharov's polylogarithms ji[ with sixth roots of unity 
as arguments. One of them is M 4j5 in Fig. 1 and the other two integrals are 
not master integrals and reduce to some master integrals in this figure. We 
performed calculations up to transcendentality weight twelve [31] and observed 
only MZV in results. 

In this paper we continue our experimental investigation of the four-loop 
massless propagator diagrams and present results for all the master integrals 
up to transcendentality weight twelve. Our motivation is twofold. First, 
we would like to check explicitly whether there are only MZV in results. 
Second, we would like to demonstrate further the power of the " dimensional- 
recurrence-and-analyticity" (DRA) method [6J| that we use. The method is 
based on the use of dimensional recurrence relations (DRR) Q and analytic 
properties of Feynman integrals as functions of the parameter of dimensional 
regularization, d, and was already successfully applied in previous calcula- 
tions SHHQHB- 

A necessary condition of the application of DRA method is the possibility 
to make an integration-by-parts (IBP) [lij reduction of integrals participat- 
ing in DRR to master integrals. To do this, we use the C++ version of the 



code FIRE [141 ]. To reveal analytic properties of the master integrals we used 



a sector decomposition [151 . 1 161 Il7j implemented in the code FIESTA [17], [18 
To fix remaining constants in the homogenous solutions of dimensional re- 
currence relations it was generally not sufficient (contrary to our previous 
work 0] on the evaluation of three non-planar diagrams) to use analytic re- 
sults for the four-loop massless propagators master integrals [l] (confirmed 
numerically by FIESTA in Ref. ]j| where one more term of the e-expansion 
was obtained) up to transcendentality weight seven. To obtain additional 
information here we applied the method of Mellin-Barnes (MB) represen- 
tation 20, [21], 22]. For each integral the number of terms calculated with 
the help of MB representation was sufficiently large to provide at least one 
consistency check of our results. 
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At the final stage of the method, we applied the PSLQ algorithm 23J as 



well as the code HPL [24| for dealing with MZV. Since the rational coefficients 
at transcendental numbers turn out to be quite cumbersome we applied PSLQ 
with the rather high accuracy of 1500 digits. 

In our results presented below, we tried to reveal as much as possible 
the homogeneous transcendentality. This is achieved by pulling out suitable 
rational functions of e and by considering a linear combination, with rational 
coefficients, of a given master integral and some of its lower master integrals. 
However, for the last three integrals, the task of finding such linear combi- 
nation is quite complicated because of the large number of the lower master 
integrals. For these three integrals we have succeeded to find their repre- 
sentation in the form of a sum of several homogeneous terms with rational 
coefficients in d. This representation allowed us to pull out a rational factor 
so that any transcendental number appeared only in the limited number of 
consecutive terms of e expansion. We choose the following loop integration 
measure 

r(d-2) d d i 



T(3 - d/2)T{d/2 - l) 2 m d l 2 

so that M3 ; i(4 — 2e) = e~ 4 . This choice corresponds to the normalization of 
Refs. 0,3.' 
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Figure 1: Master diagrams for four- loop massless propagators. Here c.l. means the com- 
plexity level — see Ref. Q. 
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+ 1^1 - 288 84 C4 C 5 + 13200 C 3Ca - + 6792 C| C 4 + ^ 

+ 495583CC, _ 1038970CO _ 882^3 + + 



+ 908472C 3 C 5 - 395297C 8 + 167994C 2>6 + 140784^ - 37260C 4 Cs + 233700C 3 C 6 

166790C 9 >9> 8380551G 2 26770505C 3 C7 108421795Cio 
+ + 109866CIC4 + — ^ + 28-^ 5^ 

- 4631 2 4 f C ^ + 390194CIC5 - 24389 f C5C6 + 488786 C4 C 7 - 
+ 3792685GC9 - 100223975(11 - 11366C 3 C 2 ,6 + 205010C 2 ,i, 8 V 



1085348C; 



+ ^1388592C 3 C5 - 324240Cs + 211248C 2)6 + 3—^ + 354744C 4 C5 

, 14nnQ o n , , M 45002569C9 M fi „, os , v 15929910(1 , 43860237C 3 C 7 
+ 1400980C 3 C6 H — h 633528C 3 C4 H Tj 1 j 

547461 15Cio 6015150C3J 2 16554412C 5 C 6 
Tj Tj + 5634364C 3 C5 - 

28467331C 4 C7 > > > > 2821800989Cn 
+ - 11209448C 3 C8 + 53532932C 2 C 9 ^ — 
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- 1436 6 8« 2 , 6 + 2893672^, - _ 4241552C , C4C5 _ 

7951810 2 256544080 254930897C 5 C 7 714442631C 3 Cg 
+ g C 3 Ce 21 ^ 2 ^ 7 + g + 27 

12873185340379Ci 2 3115920 43891225C 3 , 9 



597024 7 27 

2077280 



^2,1,1,8^ + O (e 6 ) j , 



(20) 



(l + 3e)(l + 4e) ^ 

j - 3()Sc| - 50(, () 123C 3 C4 + 514C 7 + 4862C 3 C 5 - + 1 5G(K 2 ,, ) r 



744639C 8 oooonA 1526C : 



3 



+ ^-924C 3 C 4 - 1500C 7 + 68636C 3 C 5 g— 22 + 23220C 2 ,6 + 

- 2103C4C5 + 4325(3(6 + + ( 235200C 3 C 5 710311(8 



36 J \ w " 2 



+ 85536C 2 ,6 + ^f^- ~ 36366C 4 C 5 + 55695C 3 C 6 + r^Zl^ + 2289C 3 2 C 4 
1341143CI 3816969CsC7 7815019Cio 500565^^3 , / 61040 C 3 3 



56 56 112 56 



:),.(: 



1 «nzLi c/- /- 1 1 Ri 8fin ^ ^ 460411 C 9 ^ „ n79 > 2 , ^ 60035137C 3 C 7 

- I6O4I6C4C5 + 161860C 3 C6 „ h 33072C3C4 H 

9 56 

12859479CI 134815227Ci 7724781C 3 , 7 A ro**o^ , 280574047Cn 

H : 453668UC5 H 

56 112 56 S3S 64 

+ l^Ce _ 46^. + 1309878C3Cg _ 27492nC2Cg _ 87752C3C2 6 

- 148606^)^ + (91560CIC4 + + - 

7239597C 3 , 7 ™ fi74 , 2 , ^ 3338505C 5 C 6 , 1752859 2 
^ 5183674C 3 C5 + 7, 1597650C 2 ,i,8 ^ — C 3 Ce 

, 62741559^8 ont ™ OKA A , 2990096591Cn 29828681054659Ci 2 
+ 4 29556525C 2 C 9 + ^ 2^88096 

53503881C 4 C 7 142150835Cio 170Q/1 __^ . 25832277 2 

g ^ 1 < ^ y4t > I UU^5 ^ C,2C,5 
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157544998 ^ 311533051^7 ^ 1792012205C 3 C9 ina „„ n > , 
2j C2C3C7 + ^ + 1085340C 3 C2,6 

t 1913502, t 105698899C 3 , 9 1275668 \ 5 
- 227636C 4 C2,6 H = (2(3,7 ^ g (2,1,1,8 Je 

+ 0(e 6 )J. (21) 

One may notice that rational coefficients in front of (12 involve the very 
big prime number 691 in the denominator. This is because the same num- 
ber appears in the numerator of (12 = 691tt 12 /638512875. If we use tt 12 , 
instead of (12, this number will disappear, and the highest prime factor in 
the coefficients denominators will be only 13. 

Our results are in the full agreement with the results of Ref. [l| up to 
terms considered in that paper. Partly, this is because we used some of 
the data from Ref. Q for the determination of the homo geneous parts of 
the solution of DRR. Nevertheless, our calculation can be considered as a 
nontrivial check of the results of Ref. [l| because of many constraints on the 
terms of e-expansion provided by the DRA method and fulfilled by the results 
of Ref. d. 



We also confirm the numerical results of Ref. 19 where one more term fas 



compared with Ref. [l() of the e-expansions of the most complicated thirteen 



master integrals was obtained using FIESTA [171 . 118 . 

We observe that only MZV are present in our results for the four-loop 
massless propagator master integrals. Since any other four-loop massless 
propagator integral, with any integer powers of numerators and propagators, 
can be represented, due to an IBP reduction, as a linear combination of the 
master integrals, with coefficients which are rational functions of d, we come 
to the conclusion that any four-loop massless propagator integral has only 
MZV in its epsilon expansion up to transcendentality weight twelve. This 
means that if we want to find something beyond MZV in four-loop massless 
propagator diagrams we have to go to higher transcendentality weights. This 
is certainly possible within our approach. In fact, we have chosen weight 
twelve because it looks to be already a sufficiently big number. We can 
perform calculation up to higher weights and will do our best on demand, 
for example, if somebody has reasons to believe that unusual transcendental 
constants can appear at some specific weight. 

This possibility is due to an important feature of the method that we use. 
Once we have a result in terms of a multiple series which always appears to be 
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well convergent, going to higher powers of e is an easy procedure, in contrast 
to other methodsn 




Taking our results into account one obtains more motivations to try to 
prove that there are only MZV in massless propagator diagrams. Another 
alternative is to continue to look for unusual constants in higher loops. Keep- 
ing in mind the dramatic progress of the last years in the field of evaluating 
Feynman integrals, this also looks to be a possible scenario. All results pre- 
sented here are available in computer-readable form on 
www-ttp .particle .uni-karlsruhe . de/Progdat a/ttpll/ttp 11-20/. 
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